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SEQUENTIAL PROBABILITY RATIO TEST FOR COLLISION 
AVOIDANCE MANEUVER DECISIONS 

J. Russell Carpenter* and F. Landis Markley t 


When facing a conjunction between space objects, decision makers must chose 
whether to maneuver for collision avoidance or not. We apply a well-known deci- 
sion procedure, the sequential probability ratio test, to this problem. We propose 
two approaches to the problem solution, one based on a frequentist method, and 
the other on a Bayesian method. The frequentist method does not require any 
prior knowledge concerning the conjunction, while the Bayesian method assumes 
knowledge of prior probability densities. Our results show that both methods 
achieve desired missed detection rates, but the frequentist method’s false alarm 
performance is inferior to the Bayesian method’s. 


INTRODUCTION 

We often hear it said that “space is big” by which we think decision makers mean that collisions 
between space objects are so rare that many collision avoidance maneuver (CAM) decisions are the 
result of false alarms. Nevertheless, the consequences of a missed detection might be catastrophic. 
It is our understanding of current CAM decision procedures that they do not explicitly account for 
the costs and risks of false alarms and missed detections. It seems to us that a utilitarian decision 
making procedure that balances the costs and risks of false alarms and missed detections would 
provide some benefit to decision makers. We propose herein such a procedure that is applicable to 
a common class of CAM decision problems. 

The remainder of this paper is organized as follows. In the next section, we review how the 
full CAM decision problem reduces under certain assumptions to the comparison of hypotheses 
concerning parameters of Rician densities; assumptions on which at present, our proposed method 
depends. Then, we briefly discuss collision avoidance as a decision problem involving equalizing 
utilities based on cost and risk. Next, we review Wald’s Sequential Probability Ratio Test (SPRT), 1 
which forms the basis of our proposed method. While the SPRT involves a likelihood ratio between 
two simple hypotheses, the CAM decision problem as we define it involves a comparison among 
two continuous sets of hypotheses. In the subsequent section, which we believe involves the main 
technical contribution of this paper, we discuss two methods to reduce the compound hypotheses 
to simple ones. In the penultimate section, we illustrate and compare our methods using a static 
example based on the Magnetospheric Multi-Scale (MMS) mission, and finally we conclude with a 
summary and suggestions for follow-on work. 
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PROBLEM STATEMENT 


The Simplified CAM decision problem 

Our approach to the CAM decision problem starts with the ideas of Foster and Estes. 2 In their 
study of the need for CAMs by the International Space Station, they showed that under conditions 
typical of many conjunctions, only the projection of the three-dimensional relative position error 
ellipsoid onto the conjunction plane at the time of closest approach need be considered when com- 
puting the probability of collision. Akella and Alfriend 3 later showed that a more general problem 
explicitly considering collision probability as a function of time throughout the conjunction reduces 
under certain conditions to the static problem considered by Foster and Estes. Chan 4 later showed 
that in many cases, the problem could be further simplified to that of evaluating Rician densities. 5 
In this section, we summarize the results of the cited works. 

Let { Ri , V^}, i = 1,2 represent the inertial position and velocity of two space objects. Then, the 
relative position and velocity between the objects in a frame fixed to the orbital velocity and orbit 
normal vectors of object 1 are 


r = R 2 - Ri (1) 

v = — V\ — u>i x r (2) 

where is the instantaneous orbital angular velocity of object 1. The relative position is uncertain 
due to knowledge errors, e, 

r = f + e (3) 

where f = E[r] is the nominal relative position vector. The associated error covariance matrix is 

P — E[(r — f)(r — r) T ] (4) 

= PrI + Pr2 - PR1R2 - PriR2 ( 5 ) 

where 

Pm = E [(Rt - E [Ri])(Ri - E[i^]) T ] (6) 

P RiRj = E [(Ri - E [Ri})(Rj - E [Rj}y] (7) 

The cross-covariance matrix PmRj is often negligible unless common measurements such as ranges 
between the objects have been used in the estimation process. 

Figure 1 depicts a dangerously close approach between two satellites flying in formation. The 
large red ellipsoid in the figure represents the relative position knowledge uncertainty, which is the 
propagation of post-maneuver state error uncertainty over approximately half an orbit period. One 
can see from the figure that the largest component of the knowledge uncertainty is generally along 
the velocity direction; however, there is a significant projection of the uncertainty into the conjunc- 
tion plane , i.e. the plane containing the relative position vector at the time of closest approach. The 
figure shows this projection as a blue empty ellipse in the isometric projection on the left. On the 
right, the figure shows the view normal to the conjunction plane, and the overlap between the knowl- 
edge error ellipse (“one sigma”) and the green sphere enclosing the combined hard-body radius of 
the two spacecraft. 

For such encounters, the probability of collision at the time of closest approach is the triple inte- 
gral over the hard-body volume of the density function associated with the knowledge error. The key 
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Figure 1 Stress-case simulated conjunction geometry for two MMS satellites near 
perigee. The green sphere encloses the combined tip-to-tip radius (120 m) of the two 
satellites. The red ellipsoid represents the relative position knowledge error covari- 
ance. The VNB coordinate frame on the left is aligned with the orbit velocity, orbit 
normal, and orbit bi-normal directions at the time of time of closest approach, and the 
view on the right projects the encounter into the conjunction plane. Units are meters. 

simplifying insight of Foster and Estes, 2 which was clarified and refined by Akella and Alfriend, 3 
is that the probability of collision reduces to an integration only within the conjunction plane, un- 
der certain conditions. Akella and Alfriend 3 specify the conditions under which this simplification 
holds during the short interval of the encounter, as follows: 

1 . The trajectories of the objects are straight lines. 

2. Velocity uncertainty integrates into negligible position uncertainty, and consequently the po- 
sition uncertainty is constant. 

3. The knowledge error density is Gaussian. 

Figure 2 illustrates the procedure whereby Chan 4 converted the double integral over Gaussian 
densities in the conjunction plane into a single integral of a Rician density governing the scalar miss 
distance. As Papoulis 5 describes, a Rician sample can be viewed as the root-sum-square of two 
samples drawn from Gaussian densities with the same mean and variance. Viewing the Gaussian 
samples as coordinates in a plane, the Rician sample can be viewed as a polar radius coordinate; the 
corresponding angle coordinate is free. As Figure 2 shows, Chan’s idea is to re-scale the coordinates 
of the conjunction plane so as to circularize the knowledge error ellipse. Now, the re-scaled miss 
distance has an isotropic density, and the Rician density applies. However, the region of integration, 
i.e. the hard body disk, is now elliptical, making its evaluation more difficult. To avoid this compli- 
cation, Chan suggests that we replace the now-elliptical hard body region with a circle of equivalent 
area. The result is that the probability of collision, P c , is approximated by 



(8) 
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where R represents the re-scaled hard-body radius, and I 0 denotes the modified Bessel function 
of the first kind, of order zero. We identify the Rician’s non-centrality parameter, v, with the true 
re-scaled miss distance, and view samples drawn from the Rician as representing observations of 
the re-scaled predicted miss distance. Figure 3 depicts two examples of the Rician density function. 
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Figure 2 Conversion to Rician Problem. The left pane shows the conjunction geome- 
try in the original coordinates (meters); in the middle pane, the coordinates have been 
compressed along the major axis of the knowledge error ellipse so as to circularize it, 
and the hard-body sphere has therefore been squashed into an ellipse; in the final 
pane, the elliptical hard-body has been replaced by a circle of equivalent area. 


CAM as a Decision Problem 

Assume that we receive a sequence of independent observations Z n = {zi,Z 2 , of the 

re-scaled miss distance, each with known re-scaled variance, erf. We wish to design a procedure 
that will use the observations to inform a decision concerning the need for a collision avoidance 
maneuver. In this procedure, we will compare two hypotheses, Hq and H\. The null hypothesis, 
Hq , is that the miss distance is unsafe; that is, that the true re-scaled miss distance, v, is less than or 
equal to the re-scaled hard-body radius, R. The alternative hypothesis, Hi, is that v > R. Figure 4 
illustrates the two hypotheses as families of Rician densities. 

Since we have to make this decision prior to the time of closest approach, based on a prediction 
of the miss distance, there are four possible outcomes: 

1 . A correct decision to maneuver: our miss distance prediction indicated an unsafe conjunction, 
and an unsafe conjunction really would have occurred. We will call this condition a true 
alarm , or a true positive, and denote the associated probability of detection as Pd- 

2. An incorrect decision to maneuver: our miss distance prediction indicated an unsafe conjunc- 
tion, but the miss distance was actually large enough to be safe. We will call this condition a 
false alarm, or a false positive, and denote the associated probability of false alarm as Pf a . 
This type of error creates a nuisance, but is not catastrophic. 
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Rician Densities 



Figure 3 Two examples of Rician densities, with a = 1000. For large v, the Rician 
begins to resemble a Gaussian. 



Figure 4 Families of Rician Densities for a = 500. The red Ricians correspond to 
H 0 : v < R, and the blue Ricians to Hi : v > R, where R = 50. 
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3. A correct decision not to maneuver: our predicted miss distance was large enough to be safe, 
and the actual miss distance was also large enough to be safe. We will call this condition a 
true dismissal , or a true negative , and denote the associated probability of rejection as P r . 

4. An incorrect decision not to maneuver: our predicted miss distance was large enough to be 
safe, but the actual miss distance was small enough to be unsafe. We will call this condition a 
missed detection , or & false negative , and denote the associated probability of missed detection 
as P m d • This type of error may be catastrophic. 

In principle, it should be possible to determine the cost of performing a CAM, c m , and it is also 
possible to determine the value of the objects involved in the conjunction, V . Thus, it is possible 
to associate a cost with each of the outcomes just described. Table 1 illustrates the mapping of 
decisions and conditions to outcomes and decision costs. Here, we have assumed that the cost of 
a rejection is minimal, although perhaps not zero, since there may be a cost associated with CAM 
operations regardless of whether a maneuver occurs. The expected cost of each outcome is therefore 
the product of the cost each time it occurs, and the probability with which it might occur. 


Table 1. Decision Matrix. 



Decision 

Accept Ho & maneuver 

Reject Ho & dismiss 


H 0 is true: 

True Positive: 

Missed Detection: 

8 

collision 

Probability = P d = l- P md 

Probability = P md 


will occur 

Cost — c m 

Cost = V 

a 

o 

Ho is false: 

False Alarm: 

True Negative: 

U 

trajectory 

Probability = Pf a 

Probability = P r = 1 — P/ a 


is safe 

Cost — C"ni 

Cost = e 


Since the probabilities that Ho is true, and false, are P c , and 1— P c , respectively, the total expected 
cost, c, of employing this decision procedure is 


E[c] = P c [( 1 - Pmd)Cm + P md V ] + (1 - Pc) [Pf a Cm + (1 - PfM (9) 

Neyman-Pearson theory 6 shows that for any particular fixed number of observations, there exists a 
Pareto optimal relationship between P m d and Pf a ; that is, if for example one fixes P m d = P^ d , 
then there exists a corresponding minimal P^ a (P^ d ). Letting Pj* = dPJ a /dP^ rld , optimization of 
Eq. 9 shows that the minimal cost occurs when 



Pc (V-c m ) 
(1 Pc) (Cm t) 



( 10 ) 


Since in general, one cannot arbitrarily decrease P m d without increasingP/ a , we can assume that 
P fJ P md) haS the form 

P fa(Kd) « ~k/P* md (11) 

where k is some given constant; this implies that P fa ~ P fa/ P md • Combining this with Eq. 10 
results in the approximation that 


PcPmdV ~ P/a c m(l Pc) (12) 

which shows that we minimize CAM costs by choosing P m d and P/ a such that the expected costs 
of missed detection will approximately balance the costs of false alarms. 
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WALD’S SPRT 


Let f(Z n \Ho) be the probability density associated with having obtained the set of sequential 
observations Z n given that the null hypothesis, i.e. that the miss distance is unsafe, were true*. 
Let f(Z n \Hi) be the probability density associated with obtaining the same set of observations 
given that the alternative hypothesis, i.e. that the miss distance is safe, were true. Wald’s SPRT 
compares the ratio of these probabilities to thresholds related to the false alarm and missed detection 
probabilities, as follows. 1 

Define the likelihood ratio as: 


/CZnllfr) fjZnt > g) 

" f(Zn\H 0 ) f(Z n \u <R) 

Wald recommends that we establish decision limits A and B such that whenever B < A n < A we 
will, if possible, seek another observation. If A n < B, then we will accept the null hypothesis and 
maneuver. If A n >A then we will accept the alternative hypothesis and dismiss the alarm. Wald’s 
explanations for the thresholds A and B are that we will accept the alternative hypothesis if it is A 
times more likely than the null, and accept the null hypothesis if it is 1/B times more likely than 
the alternative. Wald shows that such a procedure will terminate with probability one, and that 


A < 


p ^ a and B > 

* md 


Pfa 

1 “ Pmd 


( 14 ) 


which implies that the Pf a and P md resulting from such a procedure will be constrained by 


Pfa < B(l - P md ) ( 15 ) 

Pfa < 1 - PmdA ( 16 ) 


Note that unlike a Neyman-Pearson test which employs a fixed number of observations, the number 
of observations required for the SPRT to terminate is a random variable. 

To see how we could apply Wald’s procedure to the CAM decision problem, suppose we knew 
that the re-scaled true miss distance, v , could take on only two definite values: v o = 0 and v\ = 6R. 
Then, our hypotheses would be Hq : v — v§ — 0, and H\ : v = v\ = 6 R. Since the observations 
Z n are independent, Eq. 13 becomes 

* f{ZnW = Vl) n7 x 

" f{Zn\v = V 0 ) Ui=l^Wo,(Ti) 1 J 

from which we can see that if the z% are more consistent with v\ than with uq, the product in the 
numerator of Eq. 17 will tend to dominate the product in its denominator, and the test will terminate 
with A n > A , and vice versa. Wald points out that it may be numerically superior to use the 
log-likelihood ratio, rather than Eq. 17, 

n n 

log = T>g r ic pi I - y^logric(^|i/ 0 , CTi) (18) 

i=l i= 1 

and we employ this suggestion where applicable below. 

*Here, and subsequently, we use the function /(•) to refer to a generic probability density function. 


7 



In reality, we do not have simple hypotheses of this sort; so, in order to make use of Wald’s 
procedure, we must reduce the compound hypotheses Ho : v < R, and Hi : v > R, which Figure 4 
illustrates, to simple hypotheses. That is, we must reduce the family of red Ricians in Figure 4 to 
a single hypothesis that we can use in the denominator of Eq. 13, and reduce the family of blue 
Ricians in Figure 4 to a single hypothesis that we can use in the numerator of Eq. 13. In the next 
two subsections, we propose a frequentist and a Bayesian method to accomplish this reduction. 


Frequentist Approach 

Suppose we find the v\ E (72, oo) that maximizes the product in the numerator of Eq. 17, and also 
the vo G [0, R] that maximizes the product in the denominator of Eq. 17. Using these maximum 
likelihood estimates, \>\ and vq , respectively, in each region would accomplish our goal of reducing 
the compound hypotheses Ho : v < R and Hi : v > R to the simple hypotheses Ho : v = i>o and 
Hi : v = v\. A simple, bounded, line search will produce such estimates. 

Next, we need to determine the decision ratios A and B. Suppose we define critical measurement 
thresholds, za and zb , such that we would accept the alternative hypothesis, v > R, if a single 
measurement zi exceeds za, and we would accept the null hypothesis, v < R, if zi < zb- We 
would want to be large enough that the probability is P m d or less that the null hypothesis were 
true when zi > z a- Refer to the upper pane of Figure 5 for an illustration of this idea. Thus, 

z A = max ric _1 (l — P m dV',°l) (19) 


where 


£* = 


rz * 

= vic~ 1 (P\v, a) = f(P) : P — / ric(z\v, a) dz = 0 

Jo 


( 20 ) 


Similarly, we would want zb to be small enough that the probability is Pf a or less that the 
alternative hypothesis were true when zi < zb , as the lower pane of Figure 5 depicts. Thus, 


Zb = max ric 1 (P fa \v,ai) 

ise{R,o o) 


( 21 ) 


Then, to find A and B we merely evaluate Eq. 17 at za and zb respectively: 

A F _ max vg(fl ,oo) ric(zA|*', Cl) 

max i/6[0,fl] r ic {za\v, ci) 
b f _ max y g(fl,oo) ric(z B \v,<ri) 
max v€ [o ( jj] nc{z B \v, ci) 


( 22 ) 

(23) 


Thus, the frequentist version of our decision procedure consists of the following methodology. 
As each predicted re-scaled miss distance becomes available, accumulate the log-likelihood ratio: 

n n 

logA^ = max logric(^|^, cri) — max y^logric(^|^, cr*) (24) 

ve(R,oo ) " ve[o,R] “ 

Then, if log > log A F , dismiss the conjunction; if log A F < log B F , consider performing a 
CAM; otherwise, if log B F < log A F < log A F , attempt to obtain another observation. If it is not 
possible to obtain another observation, then it is presumably preferable to perform the CAM, with 
the understanding that the false alarm probability will be increased. 
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Figure 5. Critical measurement thresholds for frequentist Wald test. 
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Bayesian Approach 

Suppose there exists a prior density for the re-scaled true miss distance, /(V). This may be 
the case, for example, if the conjunction occurs between two formation flying spacecraft whose 
trajectories have been designed to approach to within a specified distance, z/*, which has some 
standard deviation cr*, reflecting the precision to which the formation can maneuver to achieve the 
design goal. Then, Bayes’ Rule rule says that the probability density for a single measurement being 
zi conditioned on the distance of closest approach being less than or equal to R is 


f( \v<R) = f( z = Zl ’ v ^ R ) = /tT nc(zi\v,<Ti)f(v)dv 
[1 ' ~ j Pr (v<R) f*f(, v )du 


(25) 


The probability density conditioned on the distance of closest approach being greater than R is 
similarly found to be 


f(zi\u > R) = 


f(z = Z\,v > R) 
Pr (y > R)) 


Jr ric(zi\v,ai)f(v)dv 

Jr f( u ) du 


(26) 


The integrals in Eqs. 25 and 26 must be evaluated numerically, and the infinite upper limits can be 
problematical. One way to make the integrals tractable is to assume that the prior density is also 
Rician, 

/( v) = ric(z/|z/*, cr*). (27) 

Figure 6 illustrates a particular example, and should be compared with Figure 4. 



Figure 6 Convolution of the families of densities shown in Figure 4 with a prior Rician 
density. Here, the prior re-scaled miss distance is v * = 3000, and the corresponding 
cr* = 3000, and the plot has been scaled by 1:1000 to clarify the labels. 


With this assumption, the integral in the denominator of Eq. 26 can be found using Marcum’s Q\ 
function. 7 The numerator can be expressed as the integral over all v minus the integral from 0 to 
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R, and Eq. 6.633.4 in Gradshteyn and Ryzhik 8 provides an analytic solution to the infinite integral. 
With the substitution [9, Eq. 9.6.3] Jo(iz) = Iq(z), it gives 

x exp(—ax 2 )Io((3x)Io( , yx) d x = ^ exp ^ J 0 (28) 

Straightforward algebra then leads to the intuitively pleasing result that the convolution of two 
Rician densities is a Rician density with variance equal to the sum of the variances of the component 
densities, and in particular 

ric(zi |i/, (Ji) ric(r/|z/*, cr*) d v = ric(^i|^, \J°l + (29) 

The ratio of the probabilities of Eqs. 26 and 25 then gives the likelihood ratio for a single measure- 
ment as 



f(zi\v > R) _ 

ric(zi|i/*, sjal + af) 1 

J 0 R ric M 

is*, cr*) dv 

f(zi\v < R) 

_/ 0 H ric(^i|i/, cri)ric(i/| 

\v*, cr*) dv 

J 0 °°nc^ 

|z 4 , O'*) d v 


(30) 


Assuming the measurements are independent gives the likelihood ratio for the set of n measurements 
as 


ric(^|i/*, ^J<t 2 + of) 

Jo nc(zi\y, crj) ric(i/|i/*, a*) d v 

(31) 

Finding the decision ratios A B and B B in the Bayesian approach parallels the frequentist ap- 
proach in basing these values on the results of a single measurement. We choose za to give the 
desired probability of a missed detection, 


a B _ /(gnjg ^ R) 

n ~ f{Z n \v<R) 


/ 0 ^ric(z/|z/*,er*)di/ 
J 0 °° ric(z/|z/*, a*)&v 


n 

J i= 1 


Pmd = Pr(* > z A \v < R ) 


f™ Jo nc(z\v, < 7 \) ric(/y |/y*, a*)dvdz 
f R ric(i/|z/*, o'*) d z/ 

fo A Jo 1 ric (- z l I/ ) °"i) ric(i/|i/*, o-*) d i/ d z 

J- ^ 

J 0 ric(z/|i/*,cr*)di/ 


and choose zg to give the desired probability of false alarm, 


(32) 


Pf a = Pr(,z < z B \v > -R) 


/q B JjJ ric(^|^, cri)ric(z/|i/*,cr*)dz/d,z 
J^° ric(i/|i/*,(7*)di/ 


fo B 

\v*,^/al + ol)dz- JJ B J 0 R ric(2| 

|z/, cri)ric(^ 

z/*,cr*)d^dz 

r 00 . / 

Jr ™(y\ 

|i/*, cr*) d ^ 


( 33 ) 


Then, we substitute z .4 and into Eq. 30 to find A B and B B respectively. 


11 



EXAMPLE PROBLEM 


Example Problem Description 

To test our approach, we developed a static example problem based on the MMS close approach 
described previously. For this example, the re-scaled hard-body radius, R , is 37 units, corresponding 
to 120 m for the particular close encounter geometry shown in Figure 1. The re-scaled uncertainty, 
cri, is 377 units, which corresponds to the minor axis of the unsealed 377 by 3882 meter conjunction 
plane error ellipse shown in Figure 1 . We assume that a small number of additional observations of 
the miss distance would be available subsequent to the prediction in Figure 1 . These observations 
could be more accurate, if they were based on additional tracking from assets other than those used 
for the initial observation; in any case they would be more recent as well, hence reducing error due 
to propagation. To approximate this decrease in uncertainty, we have assumed = o\ ji for this 
example. To determine the decision tresholds A and B , we used desired missed detection and false 
alarm probabilities of 1/1000 and 1/20, respectively. 

To test the frequentist approach, we generated 3000 samples of the true re-scaled miss distance 
from a uniform distribution spanning [0, 2\/3cri]. For each case, we generated a sequence of Rician 
samples, z% 9 using the true v for that case. For the prescribed missed detection and false alarm 
probabilities, the resulting thresholds were log 10 (A F ) = 2.04 and log 10 (B F ) = —.002. 

To test the Bayesisan approach, we generated 3000 samples of the true re-scaled miss distance 
from a Rician distribution assuming v * = 3000 and a * = 3000. The former is motivated by the 
unsealed design limit for close approaches for MMS of 6 km. As for the frequentist example, 
for each case, we generated a sequence of Rician samples, Zi, using the true v for that case. For 
the prescribed missed detection and false alarm probabilities, the resulting thresholds were A B = 
42.8 = 10 L63 and B b = 0.136 = 10“- 866 . 


Example Problem Results and Discussion 

We repeated the 3000-case trials several times to verify repeatable results. We determined from 
these trials that the frequentist approach reliably produces a false alarm probability of 15 d= 1%, and 
we never observed a missed detection. The average number of measurements required to resolve 
a decision, inclusive of the first observation, was two to three. In cases where the true v was very 
close to R , the frequentist approach often required many more observations; as many as 50 were 
observed. In all cases, the false alarms occurred when v « R. 

Because the frequentist approach does not assume a prior density for v, when v « R, its maxi- 
mum likelihood estimation produces estimates of v under both hypotheses that are close to v = R. 
Thus, it appears to be difficult for the frequentist approach to distinguish between the hypotheses. 
We believe this is the reason why the test fails to come closer to the desired false alarm specification 
of 1/20. Note that the limits shown by Eq. 14 apply to the actual false alarm and missed detection 
probabilities for the test, not the desired probabilities used to derive A and B, and hence they are 
trivially satisfied by the results of our trials. 

For the Bayesian approach, we also observed zero missed detections in 3000 trials, and we ob- 
served a false alarm rate of 2.8%. The average number of measurements required to resolve a 
decision was just over one, and the maximum in 3000 trials was four. 

Because the Bayesian approach has access to prior information, it appears that it is better able 
to distinguish between the two hypotheses. However, since the samples of the true v were drawn 
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from a Rician centered at 3000 units, we expect that many fewer observations near R occurred. 
Therefore, we re-ran the Bayesian case with 50,000 trials. Once again, we observed a false alarm 
rate of 2.8%, and a missed detection rate of 0. 

SUMMARY 

In this paper, we stated the case that the CAM decision problem might beneficially be viewed 
as a problem of decisions. We then applied a well-known decision procedure, Wald’s SPRT, to the 
CAM decision problem. When prior densities are available, we described a Bayesian approach to 
resolving the CAM decision problem into two simple hypotheses that applicable to the SPRT. For 
other cases, we described a frequentist approach that does not require prior information. We found 
from our examples, that, due to the nature of its reduction to Rician densities, the CAM decision 
problem is a challenging application of Wald’s theory with the frequentist approach, but that the 
Bayesian method works quite well. We believe our results show that our approaches for applying 
the SPRT to the CAM decision problem have the potential to help decision makers better understand 
the trades they face when confronted by a looming CAM decision. 

In future work, we hope to further validate and refine our method using additional MMS close ap- 
proach simulations, including if possible CAM’s to avoid other resident space objects. We also plan 
to examine the sensitivity of each of our approaches to different prior densities, to the assumption 
of independence of the measurements, and to relaxation of the some of the assumptions inherent in 
the reduction of the full CAM decision problem to the evaluation of Rician densities. 
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